Engineering the quantum transport of atomic wavefunctions over macroscopic 

distances 

A. Alberti, V. V. Ivanov, G. M. Tino and G. Ferrari 

Dipartimento di Fisica and LENS - Universitd di Firenze, CNR-INFM, 
IN FN - Sezione di Firenze, via Sansone 1, 50019 Sesto Fiorentino, Italy 



PACS numbers: 

The manipulation of matter waves represents a mile- 
stone in the history of quantum mechanics. The first 
experimental validation of matter wave behavior oc- 
curred with the observation of diffraction of matter by 
crystals pQ, and then with grating and Young's double-slit 
interference with electrons, neutron, atoms and molecules 
[2 E2 HJ |5]. More recently matter wave manipulation 
has become a building block for quantum devices such as 
quantum sensors [6] and it plays an essential role in many 
proposals for implementing quantum computers [7J |8] . In 
this letter we demonstrate coherent control of the spatial 
extent of an atomic wavefunction by reversibly stretching 
and shrinking the wavefunction over a distance of more 
than one millimeter. The remarkable experimental sim- 
plicity of the scheme could ease applications in the field 
of quantum transport and quantum computing. 

Ultracold atomic gases trapped in optical lattices 
(large and periodic ensembles of optical microtraps cre- 
ated by interfering optical laser beams) provide ideal 
tools for studying quantum transport in different regimes 
[9l [10] and quantum many body systems in periodic po- 
tentials mm USE]. 0ne of the challenges in this field 
is to coherently transfer matter waves between macro- 
scopically separated sites. This would provide a mech- 
anism to couple distant quantum bits and ultimately 
would lead to quantum information processing with cold 
atoms in optical lattices [15 j. Recently it was demon- 
strated that spatially driven lattice potentials in the pres- 
ence of a linear potential can induce a coherent dereal- 
ization of a matter wave [16 when the driving is ap- 
plied at a frequency equivalent to the Bloch frequency 
vb-, i-e., the linear potential between adjacent sites ex- 
pressed in frequency units. The derealization occurs at 
integer multiples of vb because of the resonant coupling 
between Wannier-Stark levels within the same band. The 
resonances are characterized by a sinc(27rtA^) spectral 
profile, where t is the driving time and Az/ is the detun- 
ing of the driving from the resonant frequency. The sine 
response here arises from the influence on the tunneling 
current of the relative phase <fi between the driving and 
the site-to-site quantum phase in the broadened wave- 
function. When <fi lies between and tt the wavefunction 
expands, while when it lies between tt and 2tt the wave- 
function shrinks. In particular when <fi = 2tt the wave- 
function returns to the starting point. Such a reversible 
behavior is expected provided that the evolution of the 
wavefunction is fully coherent. 

Any mechanism introducing loss of coherence would 




Figure 1: Revivals of the spatial distribution of the atoms in 
the lattice potential under strong, non-resonant driving with 
an amplitude of 10 lattice sites peak-to-peak. In the picture 
the lattice is vertically aligned and the driving is active for 
times ranging from to 3.8 seconds. The frequency detuning 
of the driving 6 is set to about 250 mHz, which results in a 
revival time of 3.8 seconds. At T=0 the in-situ RMS spatial 
width is 31 /im, while at T=3.8 s it becomes 40 /im. The color 
scale is adapted to each picture so that visibility is maintained 
with the varying atomic densities. 



in fact lead to a non-reversible broadening. However, 
in a decoherence free regime, it should be possible to 
engineer the spatial extension of the wavefunction us- 
ing the frequency offset and amplitude of the driving as 
tuning knobs. Here we experimentally demonstrate this 
new technique of matter wave manipulation by show- 
ing that coherent derealization results in an extended 
distribution corresponding to the size of the broadened 
wavefunction. This is demonstrated by observing in-situ 
the breathing of the wavefunction under non resonant 
driving conditions, and through a self interference tech- 
nique based on time-of-flight expansion. Our experimen- 
tal findings are supported by a theoretical model with 
which, using the basis of Wannier-Stark eigenstates of 
the static Hamiltonian, we can determine analytically the 
spatial wavefunction under the action of the driving. 

In order to drive (i.e., modulate) the phase of the lat- 
tice potential, we apply a sinusoidal voltage (with fre- 
quency vpzt) to the piezo-electric transducer (PZT) that 



supports the retro-reflecting mirror of the optical dipole 
standing- wave. In a first set of measurements we test the 
spatial coherence of the broadened wavefunction by ob- 
serving the periodic breathing of the atomic distribution 
while we drive the PZT with a non-zero frequency detun- 
ing Av = vpzt — vb- Fig. [I] displays the image of the 
atomic distribution under hard driving conditions. The 
frequency detuning Av is set to about 0.26 Hz such that 
revival period is about 4 seconds and the amplitude of the 
driving is set to its maximally experimentally accessible 
value, ten lattice sites peak-to-peak. The spatial profile 
is initially gaussian, corresponding to how the atoms are 
loaded into the lattice potential. The profile then evolves 
into a more complex shape as a result of the wavefunction 
broadening p~8j [19] , but then at the revival the distribu- 
tion returns to its initial profile. Starting from a size of 
31 /im it reaches an extension larger than 1.5 mm, and 
then it returns to a size of 40 /mi. In other words, the 
distribution increases its size by a factor larger than 20 
and then returns to almost to its initial value [20] . Inter- 
mediate values of the spatial extent can be obtained in a 
reproducible and reversible way by varying the amplitude 
or frequency detuning Av of the driving voltage. 

In Fig. [2] we show the time evolution of the spatial 
extent for different values of frequency detuning Av. In 
Fig. [2^i and b, with Av equal to +5 Hz and —5 Hz 
respectively, the breathing shows a revival with the ex- 
pected period of 5 Hz, and a constant visibility on a one 
second timescale, regardless of the sign of the frequency 
detuning. Fig. ^ is similar to Fig. [2^i but with Av = 0.5 
Hz. Again the atomic distribution shows a breathing at a 
frequency equal to Az/, and from the reduction of the os- 
cillation amplitude over time we can infer a e _1 damping 
time of 28 seconds. The results presented in Fig. [I] and [2] 
can not be explained classically and show a quantitative 
agreement with the analytic expression of the wavefunc- 
tion expected for the driven potential [THJ [19] [20] . This 
implies that we manipulate, and directly observe, the 
spatial wavefunction on a length scale larger than 1 mm. 

To study the coherence properties of the modified 
wavefunction we measure the interference of the wave- 
function with itself. To this end, the wavefunction is 
expanded in time-of-flight (TOF) after an adiabatic re- 
lease from the lattice potential. We record the evolution 
for the interval between to 25 ms under various driving 
conditions. First we consider the momentum distribu- 
tion under resonant driving(i.e., Av = 0). An integer 
number of sinusoidal cycles, up to 120 (equivalent to 210 
ms when vb = 574.14 Hz), is applied to the PZT such 
that the wavefunction broadens proportionally with time. 
The momentum distribution then is probed in TOF 300 
ms after the beginning of the driving, and after switching 
off the lattice potential adiabatically on a 20 /is timescale. 
Fig. [3] shows the changes in the expansion after 15 ms 
of TOF for various broadening conditions. In the ab- 
sence of a driving voltage, the thermal sample expands 
as expected yielding a gaussian profile. When we drive 
the system, we clearly observe the appearance of a non- 
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Figure 2: Dynamics and revivals of the atomic extent under 
non-resonant driving of the lattice phase. For short times 
the minimum size corresponds to the size of the non-driven 
atomic cloud, a) detuning S = vpzt — vb = +5 Hz, b) S = — 5 
Hz, c) 5 = —0.5 Hz. The driving amplitude is set to about 
1 (a and b) and 0.2 (c) lattice sites peak-to-peak. From the 
data in c) we infer a e _1 damping time of 28 s. In each graph 
the fitted line gives, within the error bars, a periodicity equal 
tocT 1 . 



Gaussian distribution. Here we can distinguish two com- 
ponents: the first is directly related to the expansion in 
absence of the modulation, while the second one has a 
spatial extent and relative weight related to the broad- 
ening. This results from the interference among the prob- 
ability amplitude originating from the different portions 
of the wavefunction as directly checked by a simulation 
of the the ballistic expansion of the broadened wavefunc- 
tion. As expected, when the free expansion due to the 
momentum dispersion is of the order of the size of the 
broadened wavefunction, the interference pattern reaches 
maximum visibility. 

There is an additional effect that, in principle, can 
complicate the observed TOF profiles. For a static lat- 
tice (i.e., one without a driving force applied to the PZT) 
the position of the interference peak exhibits periodic dy- 
namics at the Bloch frequency vb- This result, shown in 
Fig. [4] reflects the time evolution of the site-to-site phase 
differences in the wavefunction and is confirmed by simu- 
lation of the wavefunction expansion. This phenomenon 
is in direct analogy with Bloch oscillations where the 
atoms, subject simultaneously to a constant force (i.e., 
gravity in this case) and the lattice potential, sweep pe- 
riodically the first Brillouin zone [17] . In the TOF images 
this is observed as an increase of the momentum of the 
interference peack linearly in time, until it reaches the 
lattice photon recoil and changes sign. At this point the 
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Figure 3: Effect of the driving on the ballistic expansion of 
the sample. The time-of- flight is fixed to 15 ms and a reso- 
nant driving is applied during a fixed amount of cycles. Left 
side: 2D density profile. In the pictures the vertical grav- 
ity and the lattice are horizontally oriented. Right side: the 
atomic density integrated along the transverse direction of the 
lattice. The red line is the distribution expected from the bal- 
listic expansion of the analytical expression of the broadened 
wavefunction [21] . Both the colour scale on the 2D distribu- 
tion and the vertical scale on the integrated density profiles 
are the same for the various driving conditions. 
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Figure 4: Temporal evolution of the broadened wavefunction 
phase in the absence of a driving field. In a static potential 
the combined effect of the lattice and gravity induces a linear 
increase of the site-to-site phase in the broadened wavefunc- 
tion. This is observable in the ballistic expansion signal as a 
periodic evolution of the interference peak for a fixed time- 
of-flight duration of 15 ms. From top to bottom the trapping 
time is linearly increased over one Bloch period (about 1.75 
ms). 



interference peak is split in two parts. In order to sim- 
plify the treatment and neglect this additional effect, the 
pictures of Fig. [3] and the subsequent analysis are made 
by starting the TOF at point where the interference peak 
has maximum amplitude or, equivalent ly, it is centered 
in the first band before the release from the lattice. 

We have confirmed both the appearance of the inter- 
ference pattern and its Bloch like dynamics by simulating 
the evolution in the static lattice and the free expansion 
of the analytic form of the broadened wavefunction [T8] : 

+ CO 

l n W) = e- in ' 2i:VBt e i7:{n - n ' )Avt x 

n' = — co 

x J n _ n /^sin(7rAz/t)^^^ \ri) (1) 

where \n') are the Wannier-Stark eigentates of the static 
Hamiltonian (with the gravity potential and the non- 
modulated lattice) labelled by the index n' of their posi- 
tion expressed in lattice units, \n(t)) is wavefunction of 
the state \n) subject to the driving during a time t, Q 
is a coefficient accounting for the tunneling rate among 
the lattice barriers, and J n are the Bessel functions of 
the first kind. The hypothesis that the initial state cor- 
responds to the Wannier-Stark eigenstate \n) is justi- 
fied by the fact that the trapped atoms are derived from 
a thermal sample at a temperature higher than the re- 
coil temperature and they uniformly fill the first band 
(see the Methods section) . The thermal de Broglie wave- 
length is shorter than the lattice period, so when we load 



the atoms into the lattice potential, the coherence prop- 
erties among adjacent Wannier-Stark eigenstates can be 
neglected. In Fig. [3] the density profiles we derive from 
the analytical model are in remarkable agreement with 
the experimental data. It is also worth noting that the 
Bloch-like dynamics observed in Fig. [4] results only from 
the long range coherence of the broadened wavefunction. 
In the absence of driving, on the other hand, the shape of 
the distribution after ballistic expansion does not show 
any periodicity. 

As an additional check of the agreement between the 
experiment and theory we verified that the modified dis- 
tributions in the time-of-flight profiles do not result from 
a simple selection of a given momentum class. For in- 
stance, for the conditions of Fig. [3] the overall atom 
number does not change with the number of applied cy- 
cles; the atomic distribution in this case shifts to the 
central peak owing to the constructive interference be- 
tween probability amplitudes from symmetric regions of 
the wavefunction. In addition, we verified that the wave- 
function also exhibits revivals in the interference pattern 
for the TOF distribution under non-resonant driving. In 
Fig. [5] we compare the TOF profiles taken for different 
expansion times under three different conditions: with- 
out driving, at maximum expansion, and at the spatial 
revival. The driving conditions for the last condition are 
similar to those of Fig. [3] except that the frequency de- 
tuning 5 is set to about 2 Hz and the driving lasts longer 
in order to reach a maximum broadening; as in Fig. |3]with 
80 cycles. The TOF expansion at maximum broadening 
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Figure 5: Ballistic expansion under non-resonant driving. 
The periodic dynamics of the spatial extent have their coun- 
terpart in the ballistic expansion signal. From left to right: 
in-situ distribution, 10 ms and 15 ms time-of- flight. A: no 
driving. B: maximum expansion. C: first spatial revival. In 
B the driving conditions are the same as those in Fig. [3]for 80 
cycles, except here we use a 2.28 Hz frequency detuning of the 
driving field with respect to the Bloch frequency. The hori- 
zontal displacement corresponds to the free fall in the earth 
gravity field. 



is equivalent to that obtained under resonant driving ex- 
cept that now the distribution is deliberately asymmetric 
because of the different Bloch phase at the release (same 
as in fig. [4] case 5). On the other hand, the TOF ex- 
pansion at the time when the distribution shrinks to the 
initial size, is the same as that in absence of driving. 

In summary, we experimentally demonstrates a new 
method to manipulate the spatial wavefunction of cold 
strontium atoms over distances of the order of a mil- 
limeter. This manipulation preserves the quantum co- 
herence and we find that the process is fully determin- 
istic, reversible, and in quantitative agreement with an 
analytical model. The controlled stretching and shrink- 
ing of the wavefunctions may be applied for implement- 
ing new classes of quantum logic gates for q-bits such as 



neutral atoms in an optical lattice [22j [23] . It is worth 
mentioning that neutral atom optical clocks can operate 
with atoms in conditions similar to those presented in 
this work [24[ [25] . Finally, with stretched wavefunctions 
we observe a striking dependence of the distribution in 
ballistic expansion due to site-to-site phase gradients in- 
duced by external forces, regardless of the initial temper- 
ature of the sample. This approach may find applications 
in the realization of precise force sensors based on atoms 
otherwise difficult to cool. 

I. METHODS 



Experimental setup and procedure 

The experimental setup was described previously in 
[26] . About 10 6 88 Sr atoms are laser cooled to a tem- 
perature of 1 fiK and subsequently are loaded into a ver- 
tical lattice potential produed by the dipole force of a 
Xl = 532 nm laser field. The lattice is formed by retro- 
reflecting 532 nm laser light with a mirror mounted on 
top of a piezo-electric transducer (PZT). The phase of the 
optical lattice, defined to be zero at the surface of the mir- 
ror, is modulated by applying a time-dependent voltage 
to the PZT with a maximum excursion of 10 lattice sites 
peak-to-peak. The depth of the lattice potential is typ- 
ically 8 recoil energies E R (E R = h 2 /2m\ 2 L = fex 381 
nK) along the optical axis of the trapping beam and, in 
the radial direction, it decreases exponentially due to the 
gaussian spatial profile. Since initially the atoms are hot 
with respect to the recoil temperature, when they are 
loaded into the lattice potential they occupy almost uni- 
formly the first band and, with decreasing weight, the 
higher bands. Taking into account the full three dimen- 
sional lattice potential and including the gravity force, all 
the bands except the first one have a lifetime shorter than 
100 ms. We finish the atom sample preparation sequence 
by taking advantage of this property to select atoms in 
the first band, simply by holding the atoms in the static 
potential for a time of 100 ms. The atomic distribution 
is measured by absorption imaging of a resonant laser 
beam detected on a CCD camera. 
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